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Transverse-Longitudinal Coupling by Space Charge in Cyclotrons 

C. Baumgarten 1 

Paul Scherrer Institute, Switzerland^ 
(Dated: 5 January 2012) 

A method is presented that enables to compute the parameters of matched beams with space charge in 
cyclotrons with emphasis on the effect of the transverse-longitudinal coupling. Equations describing the 
transverse-longitudinal coupling and corresponding tune-shifts in first order are derived for the model of 
an azimuthally symmetric cyclotron. The eigenellipsoid of the beam is calculated and the transfer matrix 
is transformed into block-diagonal form. The influence of the slope of the phase curve on the transverse- 
longitudinal coupling is accounted for. The results are generalized and numerical procedures for the case of 
an azimuthally varying field cyclotron are presented. The algorithm is applied to the PSI Injector II and Ring 
cyclotron and the results are compared to TRANSPORT. 
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INTRODUCTION 



I. THE SIMPLIFIED CYCLOTRON MODEL 
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There is continuous interest in the understanding of 
space charge effects in isochronous cyclotrons 1-9 . The 
area is of special relevance for the conceptual designs 
of cyclotrons as energy-efficient drivers for accelerator 
driven systems (ADS) 10 '—. Nevertheless there is no self- 
consistent first order theory of matched bunched beams 
with space charge in cyclotrons up to date known to the 
author. Bertrand and Ricaud showed that it is possible 
to set up and solve linearized equations of motion for a 
simple cyclotron model 9 , but they did not determine the 
parameters of matched beams in cyclotrons and they did 
not generalize their results in order to find a numerical 
method to compute matched beams for sector-focused 
cyclotrons. Information on the matched ellipsoid is of 
special interest for modern high-performance codes like 
OPAL, that enable to simulate high intensity beams in 
the space charge dominated regim e 12 ' 13 . 

First we describe the usual simplified azimuthally sym- 
metric cyclotron model (see for instance^ or—), but in- 
clude the linearized space charge forces. The model ne- 
glects possible effects of rf-acceleration and is restricted 
to the description of a coasting beam. A modified version 
of Teng's method to decouple longitudinal and horizontal 
motion will be presented and applied to the simplified an- 
alytical mode l 16 ' 17 . Later we describe how this method 
can be applied in an iterative numerical procedure to 
determine the parameters of matched beams with space 
charge in sector-focused cyclotrons. Finally we present 
some results as computed for the cyclotrons of the PSI 
high intensity accelerator facility. 



a ' Electronic mail: christian. baumgarten@psi.ch 



The simplest relativistic cyclotron model is based on a 
magnetic field with cylindrical and mid-plane symmetrie. 
The nominal orbital frequency of an ion beam coasting 
in an isochronous cyclotron is 



1 u q R 
uj = B = B 7 , 



(1) 



where q is the charge and m the rest mass of the ions. 
The axial magnetic field is B — B(r) — Bo 7. We define a 
cyclotron length unit a = — using the orbital frequency 
uj . The particle velocity is v = ui a r so that /3 = v/c — 
r/a. Hence the relativistic factor 7 can be written as 

7= = ■ (2) 

V 1 - ( r / a ) 

The axial magnetic field in dependence of the radius r is 
then given by 



B(r) = B c 



l + e(r) 



1 r z 



(3) 



where e is a small distortion of the isochronism (e <C 
1). The radial field increase compensates the relativistic 
mass increase and the cyclotron is strictly isochronous 
for e(r) = 0. The field index n = -g 4^ is given by 



r dB 
B ~dV 



de 
dr 



so that 



7 + r 
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The phase shift per turn Arfi = 4^ is (see for instano 
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where 8 is the azimuthal angle, u c — 6 is the real or- 
bital frequency and uj^j = lu the frequency of the 
accelerating high frequency system operated at the har- 
monic number N^- The bending radius r of a particle 
with charge q in a magnetic field B is given by r = 
and hence lo c can be written as 



v mcy_8 

Tn 7 r m 7 r 



£Birl = lMM =Wo (l + £ ) 



so that 
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2nN h 



(7) 



(8) 



Hence ^ is proportional to the change of the phase shift 
per turn with radius. 

The reference orbit corresponding to the kinetic energy 
E = m c 2 (7 — 1) = Eq (7 — 1) is a circle with the radius 



r(E) 



V1 2 



1 



7 



(9) 



We will use x = r — r(E) as the horizontal coordinate 
of a specific orbit relative to the reference orbit, y as 
the axial deviation from the median plane and z as the 
longitudinal coordinate. 



II. SOME GENERAL REMARKS 

From the assumption of midplane symmetrie it follows 
that the axial coordinate y is not coupled neither to the 
horizontal nor to the longitudinal coordinate. Therefore 
it is sufficient to describe the coupling of horizontal and 
longitudinal motion. Since we assume azimuthal symme- 
trie, the matched beam ellipse must also be azimuthally 
symmetric, i.e. constant along the orbital length coordi- 
nate s. If the matched beam is constant, then the forces 
are constant - even if space charge is taken into account. 

The vector describing the (horizontal and longitudi- 
nal) displacement of an orbit at position s relative to 
the reference orbit is x = (x, x' , z, S) T where x' = 4~ 
is the horizontal direction angle with the reference orbit 



and (5 



pa 



is the momentum deviation. We use the 



formalism of linear Hamiltonian systems as for instance 
described by R. Talman 18 including the method of A. 
Wolski to compute the matched beam matrix^. The 
equations of motion (to first order) can be written as 



x'( S ) = Fx(s), 



(10) 



where F is a constant 4 x 4-matrix. The general solution 
is 



x(s) = exp (F s) x(0) = M(s) x(0) , 



(11) 



where M(s) is the transfer matrix. The matrix exponen- 
tial M(s) is 



M(s) = cxp(Fs) = 1+Fs- 



(Fs) 2 (Fs) ; 



2! 



3! 



(12) 



Let S be the skew-symmetric matrix with S 2 = — 1 



S = 



10 
-10 
1 
0-10 



(13) 



so that - since M is sympletic - the following relation 
holds: 

M T 5M = S. (14) 

If there exists an invertible transformation matrix E and 
a diagonal matrix A such that 

F = EAE-\ (15) 

then its is easy to show that the following relation holds: 

M = E exp(As)E" 1 = EAE" 1 . (16) 

E is the matrix of columnwise eigenvectors and A is the 
(diagonal) matrix of the corresponding eigenvalues of F. 
The (imaginary) eigenvalues of F are the betatron fre- 
quencies of the orbital motion. F and M have the same 
eigenvectors - and the eigenvalues of M are the expo- 
nentials of the eigenvalues of F. The eigenellipsoid oe 
defined by 

a E = Ma E M T , 



can then be written as follows^: 

a E = -EI5E -1 S , 



(17) 



(18) 



where T> is the diagonal matrix with the eigenvalues of 
<je S - which are the emittances (apart from factors ±i). 



III. THE EQUATIONS OF MOTION (EQOM) 
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The Hamiltonian is given by 

J 2 k x K x 9 

— ^ H x 
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1 2 K Z , 2 



l 2 -hx5. (19) 
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where the canonical coordinates are x and / and the mo- 
menta are given by x' and S. The EQOM (in first order) 
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where h = - is the inverse bending radius and k x — 
h 2 (1 + n) is the horizontally restoring force. K x and K z 
represent the strength of the horizontal and longitudinal 
space charge forces^: 

k 3 (1-/) Vn _ 3qIX 



X — 



K 3 
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K 3 (1-/) 
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(21) 
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The eigenvalues of F and M are: 



A 
A 
a 
b 

n 

UJ 



Diag(z f2, — i Q, i uj, —i uj) 
Diag(e l 



= K z (K x + h 2 1 2 -k x ) 
= y/a + V~a 2 — b 
= V 'a — v 'a 2 — b . 



(22) 



Note that b must be positive to give a real-valued fre- 
quency uj and hence a longitudinally focused orbit. This 
is especially important, if we allow for a (small) field er- 
ror e(r). In this case the change of the orbital frequency 
as given by 4 s can play a significant role: 



b = 



h 2 {l + n) = h 2 (l + 
h 2 (7 2 + rf) = 
K z (K x -h%) 



1) 



(23) 



if < 0, i.e. if the radial increase of the magnetic field 
is below the isochronous field increase, then the longitu- 
dinal focusing is strengthened. 

The matrix of eigenvectors E of the force matrix F is 



E 




\K zl 2 A K zl 2 A K 



The inverse matrix E 1 is given by: 



E" 



where 



2(A-B) 



-B 



i B 
Q 

A -i4 
A 



l 

1 



(24) 



(25) 



If K x < h ^ , then the longitudinal focusing frequency 
uj is imaginary and the longitudinal beam size increases 
exponentially with s. This in fact is a surprising feature 
of this type of coupling in combination with isochronism: 
Longitudinal focusing requires a strong enough horizon- 
tally defocusing space charge force. On the other hand, 

I 



A = 
B = 



h 
h 

uj 2 +K z 



(26) 



The transfer matrix M can be computed by using eq. [TBI 
and one obtains: 



M = 



A— B 



Ac-BC 
ttBS-ujAS 
AB(QS-ujs) 
c-C 



As _ B S 

Ac-BC 
AB(c-C) 
1 _ M. 



S_ 

a 
C 

AC 

s 

Bn 



t AB(C-c) 
c AB{ujS-nS) 
Be AB(Bujs- AQS) 
4-. AC-Bc 



(27) 



where 



IV. THE EIGENELLIPSOID 



C = cos (fi s) 
c = cos (uj s) 



S = sin (51s) 
s = sin (uj s) 



(28) 



The matrix T> is explicitely given by T> = 
Diag(i £i, — isi, —ie-i, i £2; 19 - The matched eigenellipsoid 
can then be calculated using eq. [T8l 
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bu 



Au 



(29) 
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The beam dimensions are given by the diagonal elements 
of the matrix representing the eigenellipsoid: 



B — A 

l 

B-A 



' Be, I Ae z \ 
. n ~ u J 
' Aex Vl+B e z u 



)• 



(30) 



where ej = and £2 = £z are identified with the hori- 
zontal and longitudinal emittance, respectively. The ax- 



ial motion can be treated separately and one finds 



h 2 v 2 y 



K„ 



(31) 



If one considers EQ.l21"ltogether with EQ.[30]and EQ.I3T] 

then it is obvious, that this result does not enable to 
start a straightforward calculation, since the beam sizes 
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depend on the space charge forces and vice versa in an 
algebraically complicated way. But with the additional 
assumption of a spherical beam, it is possible to derive a 
4th order equation for the beam size as will be shown in 
Sec. En 



V. DECOUPLING LONGITUDINAL AND TRANSVERSE 
MOTION 

If sectors are considered, then the azimuthal symme- 
trie is broken and hence the force terms in the EQOM 
and consequently the beam ellipsoid depend on the posi- 
tion s. The vertical beam dynamics can still be treated 
separately, but is has to be taken into account that the 
periodic change of the vertical beam size influences the 
space charge factors K x and K z . The design orbit usu- 
ally has to be computed numerically as described by Gor- 
don 2 ^. If this has been done, it is possible to compute the 
transfer matrix for known starting conditions, since the 
equations of motion are known and can be integrated. 
The problem is to find the correct beam dimensions for 
a given beam current and given emittances such that the 
beam is matched, i.e. such that EQ.[T7]is fulfilled. 

Tcng and Edwards described a parametrization for 
coupled motion in two and more dimensions that allows 
to find the decoupling matri x 16 ! 17 . They called their 
method symplectic rotation. We will give a brief sum- 
mary of the method and apply it to the problem of de- 
coupling longitudinal and transverse motion. 

Given two 2x2 symplectic matrices A and B that form 
a block-diagonal (i.e. decoupled) transfer matrix Tq: 



T 



A 
B 



A{in) = cos (fix) 1 
B(ju 2 ) = cos(ju 2 ) 1 




(32) 



where at, fii and are the familiar twiss-parameters. 
Then a general symplectic transfer matrix M of coupled 
motion can be written as 



M 



M n 
m N 



R Tq R 1 . 



(33) 



where M, N, m and n are 2x2 matrices and R is the 
symplectic rotation matrix. Teng suggested to write R 
in the form 



R 

R 1 



I cos ((f)) D 1 sin (4>) 
—D sin ((f) I cos (4>) 
I cos(</>) —D^ 1 sin(0) 
D sin ((f) I cos (4>) 



(34) 





Symplectic Rotation 


"Symplectic Boost" 


Tr(M-N)/2 
D t 

D 

A 
B 


cos (2 (f>) A 
sin (2 (j)f A 2 

m + S^ n T 


cosh (2tp) A 
-smh (2 V) 2 A 2 

m + S 2 n T 


A sin (2 4>) 

M-D^m tan(</>) 
N + D n tan ((/>) 


A sinh (2 lb) 

M - D~ 1 m tanh(V0 
N - Dn tanh (ip) 



TABLE I. Comparison of the method of symplectic rotation 
with the symplectic "Lorentz boost". The difference A is 
defined by A = cos (fii) — cos (^2). The matrix 5*2 is defined 

''" ! 1 and D t is defined by Dt = 2 Det(m) + 



by S 2 = 
Tr(nm). 



-1 



tain; 



^6.17 . 

M : 

N - 

m : 

n - 



D- 1 BD sin 2 ((j>) 



A cos 2 (4>) 
B cos 2 ((f)) +DAD- 1 sin 2 (4>) 
-(DA-BD) sin(0) cos(0) 
-(AD- 1 - D^B) sm((f>) cos(0) 



(35) 



The rotation angle <fi can be computed using 

\Tr(M-N) = cos (2 0) (cos (/ii) - cos (^2)) 
2 Det(m) +Tr(nm) = sin 2 (2 (f>) (cos (fii) - cos (^ j) 2 , 

(36) 

If we apply this method to the transfer matrix com- 
puted according to EQ.HH we obtain: 

\Tr(M-N) = f±^(cos(f2s)-cos(u;s)) 
2Det(m)+Tr(nm) = - (cos (Os) - cos (los)) 2 , 

(37) 

Comparison yields: 



cos (2 (f>) 
sin 2 (2 (f) 



B+A 
B-A 
4AB 



(38) 



Owe find that A > and B > and B > A 
2 so that the method fails in the case under 



From EQ. | 
since il 2 > u 

study. The method of symplectic rotation is not generally 
applicable and has to be extended. A workaround solu- 
tion was found by replacing trigonometic by hyperbolic 
functionfQ. 



R 

R 1 



/ cosh (iji) D 1 sinh (ip) 
D sinh (iji) I cosh (ijj) 

I cosh(^) — D^ 1 sinh (%ji) 
—D sinh (tfj) I cosh(ip) 



(39) 



with Det(D) = — 1, i.e. D is not symplectic - but R is 
symplectic. Tab. U compares the formulas of the symplec- 
tic rotation and of symplectic "Lorentz boost" . Applying 



where D is a symplectic 2x2 transfer matrix that de- 
scribes the structure of the coupling. Then one ob- 



1 Formally this extension is a rotation about an imaginary angle, 
similar to a Lorentz boost in Minkowski space. 
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this method to the case under study then gives: 



VI. THE ITERATION PROCESS AND EXAMPLES 



B—A 
2VAB 
B+A 



sinh (2-0) 
tanh (2 ip) 
cosh (ip) = yj 

which can be solved, 
given by: 



IAB 



B-A 



cosh (2 t/j) 
sinh (?/>) 
tanh (tp) 



_ B+A 
~ B-A 



B-A 



( 40 ) 

The matrices A and B are then 




cos (f2s) sin(f2s)/f2 
— Dsin(f2s) cos (f2s) 

cos (us) — ABuj sin (u;s) 
sin {us) /{A B u) cos(ws) 

VXb 
i/Vab o 



(41) 



The signs in B suggest a negative 02 - but this can be 
compensated by using either a negative emittance or a 
negative frequency to, since the transfer matrix is invari- 
ant against a change of the sign of u>, while ue is not. 
The transformation matrices are explicitely given by: 



R 



B — A 



A\ 
1 1/B 
A 1 
/BO 1 / 



(42) 



R 



1 















1 


-l/B 
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-A 


1 


;) 




l/B 













Note: R and R 1 are symplectic. To make the method 
complete, we give the formula to split 2x2 matrices of 
the form of A or B: 



c s/K 
-K s c 



R-2 



Rn 



COS /I 

sin \x 



i/K -i/K 

1 1 
-iK/2 1/2 
iK/2 1/2 



(43) 



The matched beam ellipsoid can be computed for given 
omittances according to EQ.[18]as soon as the diagonali- 
zation of the transfer matrix M is known. 

The fact that we have to extend the decoupling method 
of Teng and Edwards raises the question, whether the 
description of decoupling is now complete or if there are 
other cases that require further modifications or exten- 
sions. The answer to this question requires a proper two- 
dimensional extension of the Courant-Snyder theory and 
a complete survey of all possible symplectic transforma- 
tions as presented in an accompanying paper 22 . 



We have shown for the symmetric analytical exam- 
ple that a modified version of the method of Teng and 
Edwards enables to bring the transfer matrix in block- 
diagonal form according to EQ. |3U The diagonalization 
of the block-diagonal matrices is given by EQ. 031 an d 
hence the matrix of eigenvectors E and the matched 
beam ellipsoid can be constructed. 

In order to take advantage of this method for the case 
of sectored cyclotrons, the method has to be applied it- 
eratively. The goal is to compute the properties of a 
matched beam. Assumed that the beam emittances and 
the current are given as boundary conditions, one pro- 
ceeds as follows: 

1. Provide an initial guess of the beam dimensions 
cr x (s), (Ty(s) and a z (s). 

2. Compute the driving terms of the space charge 
forces K x (s), K y (s) and K z {s). 

3. Compute the one-turn transfer matrix M(s) for all 
azimuthal angles. 

4. Compute the eigenvectors of M(s) by diagonali- 
zation. 

5. Use the beam emittances to obtain the eigenellip- 
soid <je according to EQ.rjSJ 



2, if beam sizes changed significantly compared to 
the previous iteration. 

Optionally one can start the iteration with a reduced 
beam current and increase it either during the iteration 
process or use the result of the reduced beam current 
as an initial guess for higher currents. The (speed of) 
convergence strongly depends on the assumptions about 
beam current and emittance. In the following examples, 
the process converged typically after less than 20 itera- 
tions. 



A. Example: PSI Injector II 

The PSI high intensity proton accelerator facilit y 23 i 24 
(HIPA) consists of a Cockcroft- Walton pre-accelerator, a 
four-sector injector cyclotron (Injector II, 72 MeV) and 
the eight sector Ring cyclotron (590 MeV). In routine 
operation the beam current is 2.2 mA. 

The Injector II cyclotron can be modelled reasonably 
well by a hard edge approximation of the sector magnets. 
This has the advantage that the computed matched beam 
parameters can be compared to TRANSPORT 2 ^. 2 !. We 
have computed the matched beam parameters for Injec- 
tor II. The results have then be used as input parameters 
for TRANSPORT (with space charge). The results are 
shown in Fig. [T] 
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INJECTOR II, 0.842 McV, 2.2 (mA), (L^.i l .t / m2.I5.2.15J.lS) (tt mm mrad) 





FIG. 1. Transport (re-)computation of a matched beam in 
Injector II for a circulating beam of 2.2 mA at low energy. The 
boxes indicate the position of the sector magnets SM1..SM4. 
The horizontal beam width is plotted as a negative value. The 
longitudinal beam size is plotted as 2o z . 




Longitudinal Tune v z 




1 








\\\ 











^1.75 
1.5 
1.25 
1 

(1.75 

>• 01 
0.075 

0.05 

0.025 

^ 

100 200 300 400 500 600 

E (MeV) 

FIG. 2. Top: Phase shift per turn for the ring cyclotron 
computed from the measured magnetic field map. Middle: 
Horizontal and vertical betatron-tunes of the Ring-cyclotron 
as computed by EO-code (gray solid lines), the transfer ma- 
trix method over the EO (black solid lines) and with space 
charge (dashed lines). Bottom: Longitudinal betatron tune as 
induced by transverse-longitudinal coupling. The region with 
a negative slope of the phase shift is indicated by vertical lines. 
Without further field trimming, the longitudinal focusing is 
suppressed by the term —h^r- and the longitudinal beam size 
may increase. The calculation has been performed for a beam 
current of 2.2 mA and (e x ,e y ,e z ) = (1.5,2.5,0.5) 7rmmmrad. 



starting radius R eo (9 ). Besides the pure geometry of 
the orbit one also obtains a precise value for e = 1 — ^ 
for each EO^i, which then allows to compute In or- 
der to obtain the force matrix according to EQ. EUJ the 
following additional calculations have been done: 

1. Inverse bending radius h = - = B j^^ 

o p B p 



2. Local field index 

_ p dB _ B 



dB 
dr 



dB 

7d~8 



Pi 



3. Path length element As = \/r' 2 + r 2 A9. 

4. Horizontal (vertical) focusing k x = ^p- (k y 

5. h(E,9), n(E,6), j(E) and the space charge forces 
K x , K y 
matrix F 



—n ' 



K x , K y , K z are then used to compose the force 



6. The exponential series allows to compute the trans- 
fer matrix for a (short) interval [6 ... 9 + A9] : 



F 2 

F As + —As 2 + ... (44) 



M = cxp (F As) = 1 

(Typically 3 ... 4 terms are sufficient.) 



7. Multiplication of all matrices yields the one-turn 
transfer matrix M. 

8. Diagonalize M and compute eigenellipse <je- 

9. The average beam sizes are obtained from the 
eigenellipsoid and used to recompute the force co- 
efficients. 

10. The convergence is then checked and either the next 
iteration starts with step 5) or is stopped. 

Fig.[5]shows some results of the matched beam with space 
charge for the PSI ring cyclotron, the phase shift per turn 
vs. energy in the upper, the transversal tunes in the sec- 
ond and the axial tune in the last graph. The region 
where the slope of the phase shift is negative has been 
marked by vertical lines. In this region the longitudi- 
nal tune due to space charge forces apparently becomes 
imaginary. 



B. Example: PSI Ring Cyclotron 

The PSI ring cyclotron is a separated sector 
isochronous cyclotron with 8 sectors. But since the field 
of the magnets does not fall off as sharp as it does in Injec- 
tor II, a hard edge approximation does not work equally 
well. Furthermore the sectors are much more spiralled. 
For the sake of precision we used the measured field map 
B(r, 8) to compute the equilibrium orbits (EO) on a ra- 
dial grid 28 . The radius R eo {9) of the equilibrium orbit 
(EO) is then given as a function of the angle and of the 



VII. SPHERICAL SYMMETRY 

The first item in the described algorithm is to provide 
an initial guess of the beam dimensions. In the following 
we show how to obtain a possible choice of the initial 
guess for "nearly" spherical beams. If the horizontal- 
longitudinal coupling is strong enough, the beam will 
usually be approximately circular in this plane 2 ^. The 
axial size depends mainly on the vertical emittance and 
tune. 
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From EQ.[22land EQ.[26]one quickly derives: 



Analyzing EQ.[50]one finds that 



n 2 cu 2 = K z (K x + h 2 ~/ 2 -k x ) 

Q 2 + uj 2 = k x -K x - K z 

u 2 = ^/{k x -K x + K z ) 2 -Ah?^K z (45) 

1 _ 7 2 hK z 



n 



B-A 



V( fc x-^' + ^) 2 - 4,l2 7 2 K z 



With the assumption of an isochronous cyclotron (k x 
h 2 7 2 ) and spherical symmetrie, i.e. 



G 

e 
->• 



K = K. 



Oy=CTzl 



K n = K, 



K 3 7 



(46) 



one finds 



n 2 u: 2 = k 2 

n 2 + uj 2 = k x -2K 

fl 2 -uj 2 = y/k 2 -4Kk x 

1 _ K 

B ~ A hy/l-iK/k* 

The beam size a then yields using EQ. [30j 

4 _ 2 (Bui+An) 2 
_ 4e 2 



(47) 



L-4Jf ' 



so that one obtains 



4 K 3 r 2 4 e 2 r 2 
c ^ — 



3 7 



(48) 



(49) 



If one defines ctq = 



where 

x 
a 



2 -^ £ , then EQ.HHlcan be written as 



a x — 1 = . 



170 2 



g3 V2 7 r n 

3e 3/2 ^ U ' 



(50) 



(51) 



A polynom of 4th order has 4 solutions, namely the so- 
lutions of EQ. [50] are 



1 



where the abbreviations used are defined by 



2 1 / 3 Z - 



8 3 1 



> 



Z = (9 a 2 + V3V256 + 27a 2 ) 1 / 3 > 768 1 / 6 



(52) 



(53) 



The minus sign in the square root of EQ. [52] belongs to 
a pair of complex conjugate solutions and may not be 
used. Starting from low values of a, Y monotonically 
increases starting from zero so that we have to take the 
plus in front of the square root as well in order to obtain 
positive solutions also for small values of a, so that 



1 



x = 



26V3 




(54) 



lim x(a) = 1 
lim x(a) 



a 



1/3 



(55) 



so that the function x(a) — (1 + a) 1 / 3 is a reasonable 
approximation. The maximal relative deviation of about 
—0.035 appears at a ps 3/2. Therefore we suggest to use 
the approximations 



a 



Mi 
Mi 



9L - sL\ f or 

4 32 / 1U1 

1/3 for 



- a) 



< a < 5/2 
5/2 < a 



(56) 



Using the normalized emittance e — /3 7 e, the wavelength 
2 c , the cyclotron radius a = ^ = £ and 



A = — 



the relation e c 2 = -p, the values of K 3 , a and (T can 



Mo 

be written as follows: 



K, = 



3ql\ 



a 



CT = 



20 v5t eo m c 3 7 3 
3 q fig I a 

10 yf$mc(3 2 7 3 N h 
q fiQ I a y/2 7 r 
10 Vdmc/3 2 7 3 e 3 / 2 iV h 
q ^0 I /a\ 3 / 2 



(57) 




/2ae 



EQ. [55] and [57] describe the matched beam sizes for a 
spherical beam in a perfectly isochronous cyclotron. The 
real beam sizes will usually differ from these values since 
the spherical symmetrie requires that the vertical tune 
roughly equals the horizontal tune, while in most cy- 
clotrons the vertical tune is below the horizontal tune. 
Nevertheless the equations derived for the special case of 
a spherical beam can be used as starting conditions for 
the described iterative matching procedure. 



VIII. SUMMARY 

A method has been developed that allows to compute 
the parameters of a matched beam with space charge in 
cyclotrons in linear approximation for given beam current 
and known emittances. As an example, a matched beam 
in the PSI INJECTOR II cyclotron has been computed. 
The result has been used as input to TRANSPORT and 
it has been shown that the beam is matched. 

Furthermore it has been shown that the longitudinal 
focusing as induced by the space charge forces depends 
on the isochronism of the cyclotron. In case of the PSI 
Ring cyclotron, the negative slope of the phase shift per 
turn may reduce or even destroy the longitudinal focusing 
effect. The results are important for the design of high 
intensity cyclotrons. If beam emittance and current are 
chosen accordingly, then the space charge induced longi- 
tudinal focusing allows to replace flat-top cavities by ac- 
celerating cavities and thus increase the energy gain per 



8 



turn and the turn separation significantly. This is im- 
portant to minimize beam losses and activation of com- 
ponents. The flat-top cavities of the PSI Injector II cy- 
clotron are already operating as accelerating cavities and 
a complete replacement is planne d 29 ' 30 . 
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